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AUTHORS'  SUMMARY 

This  study  is  concerned  with  an  experimental  investigation  of  the  different  ways  in 
which  turbulence  appears  in  an  incompressible  laminar  boundary  layer  on  a  swept  wing. 
Transition  is  detected  by  the  use  of  hot  films  attached  to  the  model.  The  results  of  two 
series  of  measurements  are  given:  in  the  first,  the  sensors  are  placed  along  the  chord  and  the 
recorded  signals  allow  the  analysis  of  problems  in  streamwise  and  crossflow  instability,  and 
leading  edge  contamination.  In  the  second  series  of  experiments,  hot  films  are  distributed 
along  the  span,  close  to  the  leading  edge.  Leading  edge  contamination  is  studied  in  detail 
and  possibilities  of  relaminarisation  are  considered. 
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1  INTRODUCTION 

A  swept  wing  of  infinite  span  represents  the  simplest  case  of  three-dimensional 
flow.  As  shown  in  Fig  1,  two  systems  of  co-ordinates  are  generally  defined:  one  (X,  Z,  y) 
is  related  to  the  wing,  while  the  other  (x,  z,  y)  is  related  to  the  external  flow.  In  both  cases, 
the  y  direction  is  normal  to  the  wall  and  the  other  co-ordinates  are  measured  along  the 
aerofoil,  starting  at  the  geometrical  leading  edge  and  perpendicular  to  it.  The  components  of 
mean  velocity  are  denoted  (U,  W,  v)  in  the  (X,  Z,  y)  system,  and  (u,  w,  v)  in  the  (x,  z,  y) 
system.  The  velocity  profiles  u(y)  and  w(y)  are  referred  to  respectively  as  the  streamwise 
and  crossflow  profiles.  Outside  the  boundary  layer,  we  =  0  . 


Using  this  geometry,  the  fundamental  hypothesis  is  that  derivatives  in  the  direction 
Z  of  the  span  are  zero.  It  follows  that  the  component  We  of  the  external  velocity  is 
constant  over  the  whole  wing  and  equal  to  Woo  =  Qoo  sin  tp  (see  notation  in  Fig  1).  If  the 
distribution  ue(X)  is  known,  for  example  by  measurement  of  static  pressure,  it  is  possible 
to  deduce  the  distribution  Ue(X)  from  the  equation: 

U*(X)  =  u^(X)  -  (Q^  sin  <p)2  (1) 

In  a  perfect  fluid,  it  may  be  shown  that  the  profile  Ue(X)  is  that  which  would  occur 
on  an  unswept  wing  situated  in  a  flow  Uoo  =  Q«>  cos  cp  .  In  particular,  Ue  =  0  along  the 
attachment  line,  the  distinct  flow  line  which  separates  the  flows  of  the  upper  and  lower 
surfaces. 

Even  if  a  swept  wing  of  infinite  span  has  simple  geometry,  the  process  of  transition 
may  turn  out  to  be  considerably  more  complex  than  for  two-dimensional  configurations.  In 
fact  turbulence  may  arise  from  three  different  mechanisms:  streamwise  instability,  crossflow 
instability,  and  leading  edge  contamination  (v.  for  example,  Poll1,  Amal^,  Saric  and 
Reed3). 
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Streamwise  instability  is  related  to  the  properties  of  the  streamwise  profiles  u(y) . 
Since  these  resemble  the  classical  Falkner-Skan  profiles,  it  is  similar  to  that  associated  with 
two-dimensional  flow  and  is  the  cause  of  transition  of  decelerated  flow.  However,  the 
inevitable  presence  of  a  point  of  inflection  in  the  crossflow  profiles  w(y)  may  render  these 
very  unstable  in  those  regions  where  they  develop  most  rapidly,  ie  in  the  region  of  the 
leading  edge.  Crossflow  transition  is  encountered  mostly  in  accelerating  flows. 

From  a  practical  point  of  view,  the  prediction  of  streamwise  and  crossflow  transition 
may  be  made  by  applying  simple  criteria4  or  by  the  use  of  the  linearised  theory  of  laminar 
instability.  The  latter  technique  consists  of  calculating  the  growth  of  unstable  frequencies 
and  placing  the  theoretical  transition  at  the  point  where  the  most  dangerous  frequency  has 
grown  by  a  ratio  en  ,  with  n  in  the  approximate  range  7  to  9. 

The  third  transition  mechanism  referred  to  above,  leading  edge  contamination,  is  of  a 
very  different  nature.  Taking  the  simple  example  of  a  swept  wing,  if  this  is  influenced  by  a 
solid  surface  (fuselage,  tunnel  wall),  the  turbulence  convected  on  to  this  wall  may  become 
propagated  along  the  attachment  line  and  possibly  render  the  wing  completely  turbulent. 
This  basically  non-linear  phenomenon  is  an  example  of  what  Morkovin  has  called  a 
"bypass”5,  in  the  sense  that  it  cannot  be  treated  by  the  laminar  theory  of  instability.  Under 
incompressible  conditions,  experiment  has  shown  that  the  appearance  of  contamination  is 
determined  by  the  value  R  of  the  Reynolds  number  calculated  for  the  attachment  line 

W  rj 

R  =  —  (2) 

v 

v  and  t]  represent  respectively  the  kinematic  viscosity  and  a  characteristic  length  defined 
by 

ti  =  (v/k)’/2  ,  with  k  =  (dUe/dX)x_x  ,  (3) 

p 

where  Xp  is  the  ordinate  of  the  attachment  line.  Generally,  k  is  related  to  the  local  radius 
of  curvature  of  the  wing.  It  is  a  result  of  a  large  number  of  experiments  that,  if  R  is  less 
than  245,  the  turbulence  arising  from  the  wall  is  suppressed  and  disappears  more  or  less 
quickly  along  the  attachment  line.  If  this  is  not  the  case,  there  is  contamination,  an  increase 
in  the  magnitude  of  turbulent  phenomena  which  end  up  occupying  the  whole  of  the  leading 
edge.  This  simple  criterion,  used  by  Pfenninger  since  19656,  has  since  that  time  been 
largely  confirmed  (Gaster7,  Cumpsty-Head8,  Poll9,  for  example). 

This  article  outlines  two  series  of  experiments  in  the  FI  and  F2  wind  tunnels  of  the 
ONERA  Fauga-Mauzac  centre,  carried  out  in  order  to  check  the  methods  of  calculation  of 
transition  on  swept  wings.  During  the  first  series  of  measurements  carried  out  on  a  wing 
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with  a  cambered  leading  edge  (F2  wind  tunnel)  the  three  forms  of  transition  were  studied. 
During  the  second  series  of  experiments,  conducted  on  another  wing  in  the  FI  wind 
tunnel,  an  analysis  was  made  in  greater  detail  of  the  problems  of  leading  edge  contamination 
and  of  possible  relaminarisation. 


2  EXPERIMENTS  ON  A  WING  WITH  CAMBERED  LEADING  EDGE 

(WING  No.l) 

2.1  Experimental  arrangement  and  method  of  measurement 

These  experiments  were  carried  out  in  the  F2  wind  tunnel,  which  has  a  test  section 
of  1.4  m  x  1.8  m  and  a  length  of  5  m.  A  detailed  description  is  given  at  Ref  10.  The  air 
flow  is  provided  by  a  fixed-bladed  fan  of  variable  speed,  giving  an  air  velocity  in  the  range 
from  zero  to  100  ms'1  within  the  section.  External  turbulence  is  less  than  0.1%. 

The  model  is  a  wing  having  an  ONERA  D  aerofoil  and  a  C  chord  of  0.3  m  normal  to 

the  leading  edge.  It  has  a  leading  edge  cambered  between  0  and  20%  chord.  The  presence 

of  this  leading  edge  obviously  contributes  to  pressure  distributions  which  differ  from  those 

observed  on  a  classic  ONERA  D  aerofoil.  The  wing  is  attached  to  a  half-body,  which  is  in 

turn  attached  to  one  of  the  vertical  walls  of  the  section.  Fig  2  shows  the  aerofoil  and  a 

diagram  of  the  experimental  arrangement  in  the  section. 

•* 

Three  kinds  of  measurement  were  carried  out:  pressure  distribution  at  the  wall,  flow 
visualisation  by  sublimation,  and  analysis  of  the  signals  provided  by  the  hot  films.  The  hot 
films  were  attached  to  the  model  at  between  2.5%  and  86%  chord,  and  their  signals  were 
recorded  for  a  large  number  of  combinations  of  the  three  parameters,  (p  ,  a  ,  and  Qoo  , 
respectively  sweep  angle,  geometrical  angle  of  incidence  and  nominal  flow  velocity  in  the 
section.  The  positions  of  the  hot  films  are  shown  in  Fig  2b;  they  were  chosen  in  order  to 
avoid  interactions  between  neighbouring  sensors.  The  absence  of  interference  between 
sensors  was  confirmed  by  comparing  the  transition  positions  indicated  by  the  hot  films  and 
by  flow  visualisation  at  the  wall,  and  the  results  were  shown  to  be  identical. 


x  hot  films 

Working  section  axis - 


a)  Cross-section  of  the  aerofoil 


b)  Arrangement  in  the  working  section 


Fig  2  Experimental  arrangement  in  the  F2  tunnel 
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2.2  Configurations  studied 

Three  sweep  angles  (9  =  49,  55  and  61°)  and  four  geometrical  angles  of  incidence 
(a  =  0,  -2,  -4  and  -8°)  were  studied.  Typical  distributions  of  the  external  velocity  are 
given  in  Fig  3.  On  account  of  the  cambered  leading  edge,  the  curves  show  complex 
variation  in  the  region  of  negative  pressure  gradient. 


2.3  Results  of  streamwise  and  crossflow  instability 

Analysis  of  the  results  has  shown  that,  in  a  large  number  of  cases,  transition 
occurred  as  a  consequence  of  crossflow  instability. 

For  example,  Fig  4  shows  the  variation  of  the  transition  abscissa  Xy  as  a  function 
of  the  velocity,  Qx  ,  for  <p  =  49°  ,  a  =  -2°  .  The  experimental  data  were  compared  with 
predictions  derived  from  the  separate  application  of  criteria  for  transition  developed  at 
CERT/DERAT4.  The  criterion  for  streamwise  transition  is  an  extension  of  that  developed 
by  Granville;  it  is  expressed  by  an  equation  of  the  form 


Re  Re  f(A2,  Tu)  .  (4) 

T  cr 

Rq  is  the  Reynolds  number  based  on  the  momentum  thickness  of  the  streamwise 
profile;  the  suffixes  T  and  cr  denote  the  transition  point  and  the  starting  point  of 
instability.  A2  is  a  parameter  of  the  streamwise  pressure  gradient  (mean  Pohlhausen 
parameter),  and  Tu  is  the  level  of  external  turbulence. 
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- Laminar  calculations 

O  Experiment 


i _ i _ i _ i _ i _ i _ i _ 

0  40  80  CUms-')  120 


Fig  4  Variation  of  the  transition  abscissa  as  a  function  of  CL. 

(<p  =  49°,  a  =  2°) 

Comparison  with  the  results  of  the  criteria  for  transition 

The  crossflow  criterion  is  an  empirical  correlation  between  two  complete  parameters 
considered  as  the  transition  point:  the  crossflow  Reynolds  number 


and  the  streamwise  shape  parameter  H-. 

In  the  region  where  the  flow  is  strongly  accelerated  (X/c  <  0.25),  the  results  of 
applying  the  crossflow  criterion  are  in  good  agreement  with  the  measurements.  However, 
further  downstream  where  the  flow  is  decelerated,  the  hot  films  detected  bursts  of 
turbulence  not  predicted  by  the  criteria.  As  shown  in  Fig  4,  the  transition  points  derived 
from  the  streamwise  criterion  are  situated  well  downstream  of  the  experimental  results. 
Discussion  of  this  problem  may  be  considered  from  two  aspects: 

(a)  the  criteria  were  established  from  experiments  in  which  the  transition  was  the 
result  of  instabilities  which  were  wholly  either  streamwise  (x  direction)  or  crossflow 
(z  direction).  It  is  possible  that,  between  25%  and  90%  chord,  the  most  unstable 
directions  change  progressively  from  the  z  direction  to  the  x  direction.  Only 
rigorous  calculations  of  stability,  associated  with  the  en  method,  may  then  predict 
transitions  of  an  intermediate  type.  Such  an  approach  was  adopted  by 
Cebeci  et  a/11,  and  the  results  which  they  achieved  are  to  be  presented  in  another 
paper  in  this  congress. 

(b)  The  classical  theory  of  linear  instability,  as  well  as  the  semi-empirical 
criteria,  assume  that  there  is  no  interaction  between  streamwise  and  crossflow 
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instabilities.  Theoretical  research  such  as  that  conducted  by  H.  Reed12  has  shown 
that  stationary  waves  resulting  from  crossflow  instability  may  double  the 
magnification  coefficients  of  streamwise  instabilities,  and  consequently  lead  to 
premature  transition.  These  problems  of  interaction  are  of  course  very  difficult  to 
take  into  account  in  stability  calculations.  They  may  be  introduced  in  an  empirical 
manner  into  the  streamwise  criterion  by  modifying,  for  example,  the  pressure 
gradient  parameter  by  a  function  of  R§2  .  Such  work  has  been  carried  out  at 
CERT/DERAT  and  has  led  to  an  improvement  in  prediction  *3.  However,  the  limited 
number  of  measurements  available  to  support  this  modification  does  not  allow  a 
definite  conclusion  to  be  reached. 

2.4  Results  of  leading  edge  contamination 

In  this  series  of  experiments,  it  was  assumed  that  leading  edge  contamination 
appeared  when  turbulent  fluctuations  were  detected  on  the  upstream  hot  film,  placed 
throughout  at  X/c  =  0.025  ,  although  the  abscissa  Xp  of  the  attachment  line  varied 
appreciably  as  a  function  of  incidence. 

Fig  5  shows  some  examples  of  fluctuating  signals  obtained  for  cp  =  55°  ,  a  =  -8° 
and  two  values  of  velocity  Qoo .  For  each  value  of  the  latter,  the  signals  obtained  from  the 
ten  hot  films  were  simultaneously  recorded.  Measurements  of  wall  pressure  indicate  that  the 
attachment  line  is  situated  at  X/c  =  0.05  . 


(arbitrary  scale) 

x/c 

0.025 

0.15 

0.20 

S^JL  XJL  aA 

0.25 

\a_  yv_A. 

0.30 

i  i  i 

\A_i*Nrv  rv  a.  ajv 

0.35 

0.45 

0.55 

0.67 

0.86 

0  0.04  0.08  0.12  t(s)  0.16 

a)  0*  =  33  ms-1 


X'  (arbitrary  scale)  x/c 

0.025 

0.15 

0.20 

0.25 

0.30 

0.35 

0.45 

0.55 

0.67 

0.86 

0  0.04  0.08  0.12  t(s)  0.16 

b)  Qm  z  35  ms'1 


Fig  5  Examples  of  fluctuating  hot  film  signals 
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At  the  lower  speed  (Q  =  33  ms-1,  Fig  5a)  the  region  of  the  leading  edge  remains 
laminar,  but  a  large  part  of  the  wing  is  "contaminated"  by  turbulent  spots  arising  from  the 
area  of  the  junction  between  the  wing  and  fuselage.  An  explanatory  diagram  of  the  situation 
is  given  in  Fig  6:  turbulent  structures  are  created  at  the  root  by  the  turbulent  boundary  layer 
of  the  fuselage,  and  as  R  is  too  small,  they  cannot  develop  along  the  attachment  line,  but 
those  that  are  convected  along  the  line  of  external  flow  are  detected  by  the  hot  films 
downstream.  It  is  also  to  be  noted  that  the  recordings  allow  the  signatures  of  distinctive 
spots  to  be  followed  in  the  X  direction.  However,  as  the  external  streamline  was  almost 
perpendicular  to  the  row  of  hot  films,  it  was  not  possible  to  calculate  the  convection  velocity 
of  the  spots.  It  is  noted  on  the  recordings  that,  paradoxically,  the  beginning  of  a  spot  is 
often  detected  earlier  on  film  n  than  on  film  n  -  1  upstream.  The  difference  between  the 
direction  of  displacement  of  the  spots  and  the  orientation  of  the  row  of  hot  films  allows  this 
phenomenon  to  be  explained;  for  example,  in  Fig  6  the  spot  in  state  2  will  be  detected  by  hot 
film  9  before  it  is  detected  in  state  3  by  hot  film  8. 

When  the  velocity  in  the  section  was  increased  from  33  to  35  ms'1,  the  number  of 
turbulent  spo^s  grew  rapidly  and  some  turbulent  bursts  of  air  appeared  at  the  first  position. 
Under  these  conditions,  the  Reynolds  number  R  was  equal  to  276,  which  is  taken  to  be  the 
value  for  the  onset  of  contamination. 


It  was  possible  to  observe  leading  edge  contamination  for  five  other  (tp,  a)  pairs. 
The  first  turbulent  bursts  were  detected  on  the  upstream  sensor  at  RD=  258  ±  18  and  the 
attachment  line  was  completely  turbulent  at  Rp  =  309  ±  12  .  These  values  are  in  good 
agreement  with  those  found  during  earlier  work. 


A  Spots  coming  from  the  leading  edge 
•  'natural'  transition 


=  55s.  Of  =  -8s 
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Fig  6  Interpretation  of  the 
hot  film  recordings 


Fig  7  Transition  abscissa  for  two  cases 
of  leading  edge  contamination 
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In  actual  fact,  contamination  may  appear  in  a  manner  which  differs  according  to  the 
case  under  consideration.  In  Fig  7,  the  transition  position  is  drawn  for  two  experimental 
configurations,  as  a  function  of  the  velocity  Q<„  .  (The  transition  position  is  taken  at  the 
streamwise  abscissa  at  which  the  first  turbulent  fluctuations  are  recorded).  With  9  =  55° 
and  a  =  -8°  (the  situation  corresponding  to  the  signals  shown  in  Fig  5),  Xj  rose  very 
rapidly  from  the  trailing  to  the  leading  edge  with  neither  streamwise  nor  crossflow 
instability.  The  situation  was  more  complex  with  <p  =  61°,  a  =  -2°.  When  Q*,  increased 
from  40  to  65  ms**,  the  transition  shifted  slowly  from  30%  to  20%  chord.  It  was  initiated 
by  crossflow  instability,  as  shown  by  the  application  of  the  crossflow  criterion.  At  higher 
speeds  R  exceeds  245  and  leading  edge  contamination  is  evident.  This  leads  to  a  complex 
situation  in  which  two  mechanisms  for  producing  turbulence  are  simultaneously  present. 
An  example  of  this  situation  is  shown  in  Fig  8,  in  which  the  velocity  Qoo  vas  slightly  less 
than  that  at  which  contamination  apperrs.  Some  turbulent  spots  were  produced  at 
X/c  =0.15  ;  they  occurred  near  to  the  attachment  line  in  the  region  of  the  junction  of  wing 
and  fuselage.  At  following  positions  a  mixture  of  spots  arising  from  the  attachment  line, 
together  with  those  which  are  the  result  of  a  transition  by  crossflow  instability,  may  be 
observed. 


In  the  course  of  these  experiments,  development  of  turbulence  in  the  X  direction  is 
of  particular  interest.  In  a  second  series  of  measurements  described  in  the  next  paragraph, 
attention  was  directed  towards  the  development  of  spots  along  the  span,  in  areas  adjacent  to 
the  leading  edge. 


x/c 
0025 

0  15 

0  20 

0.25 

0  30 

035 

0.45 

0.55 

0.67 

0  86 

C  0.04  0.08  0 12  t(s)  0.16 


X'  (arbitrary  scale) 


Cas 

9° 

a° 

Pt  (bars) 

A 

40 

2.5 

3 

B 

40 

2.5 

2 

C 

30 

2.82 

3 

D 

30 

11.3 

3 

E 

30 

11.3 

2 

F 

30 

11.3 

1 

G 

40 

10 

2 

H 

40 

10 

1 

Table  1  -  configurations  studied 


Fig  8  Transition  signals  for  9  =  61°  ,  Table  1  -  Configurations  studied 

a  =  -2°  ,  Q„o  =  66  ms-1 
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3  EXPERIMENTS  ON  WING  No.2 
3.1  Experimental  arrangement 

A  second  series  of  experiments  was  carried  out  in  the  FI  wind  tunnel  of  the 
Fauga-Mausac  Centre  on  the  so-called  "hyper- lift  variable  sweep"  wing  of  the  ONERA 
Aero  Department.  An  aerofoil  section  and  a  diagram  of  the  arrangement  of  the  model  in  the 
wind  tunnel  are  shown  in  Fig  9. 

The  working  section  is  10  m  long.  At  its  exit,  the  height  is  3.5  m  and  the  width, 
4.5  m.  The  tunnel  is  driven  by  a  constant- speed  axial  fan  of  variable  pitch,  so  that  the 
velocity  in  the  section  can  be  m-'de  to  vary  from  zero  to  about  100  ms-1.  The  installation 
operates  at  ambient  temperature,  but  the  total  pressure  may  be  controlled  between  1  and 
3  bars  (v.  f!4J). 


C=  0.5  m 


a)  Section  of  aerofoil 


b)  Arrangement  in  the  working  section 

Fig  9  Experimental  arrangement  in  the  FI  wind  tunnel 

The  model  essentially  comprised  a  wing  made  to  the  RA16  SCI  aerofoil,  except 
between  X/c  =  0  and  X/c  =  0.20  ,  where  an  AMD-BA  leading  edge  was  employed.  It  was 
fitted  wi:h  a  trailing  edge  flap  set  at  10°  in  all  the  experiments  (Fig  9a).  The  chord  had  a 
constant  value  of  0.5  m.  With  an  angle  of  sweep  9  =  40°  ,  the  upper  edge  of  the  wing  was 
situated  2  m  above  the  floor.  The  model  was  fitted  with  eight  rows  of  static  pressure 
sensors  along  the  span;  only  the  rows  of  sections  2,  4,  5,  6  and  7  were  employed.  Their 
location  is  shown  in  Fig  9b;  at  9  =  40°  ,  these  rows  are  parallel  to  the  floor  of  t^e  working 
section. 
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During  the  commissioning  of  the  tunnel,  measurements  of  the  boundary  layer  were 
carried  out  at  a  number  of  points  in  the  working  section.  At  the  fixing  point  of  the  model  on 
the  floor,  measurements  were  obtained  of  classical  profiles  of  a  turbulent  boundary  layer, 
with  a  physical  thickness  of  around  10  cm  in  the  range  of  Reynolds  numbers  encountered  in 
this  study.  Along  the  axis  of  the  working  section,  the  degree  of  external  turbulence  is  of  the 
order  of  0. 1  %. 

3.2  Methods  of  measurement  -  configurations  studied 

As  in  the  preceding  study,  the  hot  films  were  the  essential  means  of  detection  of 
transition.  The  transient  signals  simultaneously  provided  by  eight  sensors  were  recorded  by 
analogue  methods,  then  digitised  and  processed  at  the  CERT  computation  centre. 

Eight  configurations,  numbered  A  to  H,  were  studied.  Table  1  shows,  in  each  case, 
the  values  of  the  angle  of  sweep  (p  ,  the  geometrical  angle  of  incidence  a  and  the  total 
pressure  Pi  .  The  configurations  A,  B,  C  correspond  to  a  normal  angle  of  incidence 
an  =  a/cos  <p  of  3.26°  ,  (referred  to  as  "low  incidence"),  while  the  five  other  configurations 
correspond  to  an  =  13.05°  ,  (referred  to  as  "high  incidence"). 


line 

b)  Case  H 


Fig  10  Examples  of  pressure  distributions 
3.3  Pressure  distribution  -  attachment  line 

Fig  10  shows  the  distributions  of  the  pressure  coefficients  Cp  for  two  typical 
examples  of  low  (case  B)  and  high  (case  H)  incidence.  The  measurements  obtained  at 
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sections  2, 4, 6  and  7  are  shown  as  a  function  of  X/c  ,  where  X  is  the  curvilinear  abscissa 
measured  along  the  aerofoil  from  the  geometrical  leading  edge.  By  convention,  the  positive 
values  of  X  are  associated  with  the  upper  surface  and  the  negative,  with  the  lower  surface. 
The  increase  at  the  point  of  over-velocity  on  the  upper  surface  as  a  increases  is  particularly 
spectacular.  The  point  at  which  Cp  achieves  a  maximum,  theoretically  equal  to  cos2<p  for 
an  infinite  swept  wing,  identifies  the  attachment  line. 

By  measuring  Cp  ,  it  was  possible  to  calculate  the  variation  of  the  velocity  Ue(X) , 
and  then,  by  using  equation  (1),  the  distribution  Ue(X)  ,  the  projection  of  the  external 
velocity  perpendicular  to  the  leading  edge.  The  variations  of  Ue/Q<»  as  a  function  of  X/c 
are  shown  on  Fig  1 1  for  the  two  cases  under  consideration  in  the  previous  figure.  Only  the 
measurements  in  the  region  of  the  leading  edge  are  plotted,  and  the  velocities  on  the  lower 
surface  are  shown  with  a  negative  sign  in  order  to  enable  the  attachment  line  (Ue  =  0)  to  be 
better  visualised. 


ue/Q<o 


Section 

•  2 

00®$  * 

§xx..- 
x  •  • 

x  4 

• 

a  5 

a 

O  6 

X 

a  7 

a* 

X 

Lower 

surface 

S* 

$ 


•  |\ 


Upper 

surface 


Attachment  line 


-0.1 


x/c 


0.1 


-0.2  -0.1  x/c  0 


a)  Case  B 


b)  Case  H 


Fig  11  Distribution  of  velocity  normal  to  the  leading  edge  In  the 
region  of  the  attachment  line 

It  appears  that  the  results  obtained  in  the  different  sections  do  not  fit  entirely  onto  a 
single  curve  and,  because  of  this,  the  abscissa  Xp  of  the  attachment  line  varies  little  along 
the  span. 

3.4  First  series  of  hot  film  measurements 

As  shown  in  Fig  12,  for  cases  G  and  H,  eight  hot  films  were  first  attached  along  the 
leading  edge  of  the  wing,  at  approximately  2%  chord  on  both  sides  of  the  attachment  line  as 
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determined  by  pressure  measurements.  The  films  numbered  2,  4,  6  and  8  were  situated  on 
the  upper  surface,  while  those  numbered  1,  3,  5  and  7  were  situated  on  the  lower  surface. 
In  Fig  12  the  position  of  the  sensors  along  the  span  is  shown  in  relation  to  the  lines  of  static 
pressure  measurements. 


Fig  12  Positions  of  hot  films  and  pressure  measurement  lines 

During  adhesion  the  hot  films  were  aligned  so  that  the  sensitive  element  would  be  nearly 
perpendicular  to  the  external  flow  streamline. 

For  each  of  the  eight  configurations  under  consideration,  a  progressive  sweep  in 
velocity  was  undertaken  to  define  precisely  the  onset  of  contamination.  For  example, 
Fig  13  shows  signals  recorded  in  case  H  for  Q*,  =  35,  61  and  64  ms'1  .  At  the  lower 
velocity  (Fig  13a)  only  film  1  indicates  turbulent  fluctuations  which  arise  from  the  turbulent 
boundary  layer  in  which  the  hot  film  is  situated,  at  the  floor.  At  Q,*  =  61  ms'1  ,  a  number 
of  "spot"  type  turbulent  structures  are  recorded  by  film  2;  they  are  seen  to  diminish  in 
number  at  films  3  and  4,  then  to  remain  constant.  At  films  4  to  8  these  spots  are  seen  to  be 
enlarged;  from  then  on  they  develop  along  the  attachment  line,  indicating  the  start  of  leading 
edge  contamination.  (Attention  must  be  paid  to  the  interpretation  of  these  signals;  for 
example,  on  film  8  two  more  spots  may  be  counted  than  on  film  6.  In  fact,  the  two  spots  at 
the  left  on  film  8  also  appear  on  film  6,  but  earlier  in  the  recording).  An  increase  of  5%  in 
the  velocity  Q«,  (from  61  to  64  ms-1)  leads  to  a  very  rapid  growth  in  the  number  of  spots, 
as  a  comparison  of  Figs  13b  and  13c  shows.  At  95  ms-1,  all  signals  are  fully  turbulent. 

Next,  the  calculation  of  R  at  the  onset  of  contamination  was  carried  out  using 
equations  (2)  and  (3).  The  critical  part  consisted  of  estimating  the  parameter  k  =  (dUe/dX) 
at  X  =  Xp  .  Fig  1 1  illustrates  the  difficulty  of  obtaining  an  accurate  value  of  the  velocity 
gradient  in  the  region  of  the  attachment  line,  with  this  gradient  varying  somewhat  along  the 
span.  After  smoothing  of  the  experimental  points,  a  mean  value  for  k  was  taken  on  the 
five  sections  on  which  the  pressure  was  measured.  The  onset  of  contamination  was  located 
at  1^=  251  ±11,  the  leading  edge  being  completely  turbulent  at  Rp  =  318  ±  22  .  These 
results  are  very  close  to  those  obtained  in  the  first  series  of  experiments. 
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Fig  13  Signals  obtained  from  hot  films, 
(p  =  40° ,  a  =  10° ,  P|  =  1  bar 
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In  the  majority  of  recordings  in  which  RD  <  R  <  Rp ,  it  was  possible  to  follow  the 
movement  of  the  indivudal  spots  from  one  hot  film  to  another,  and  to  derive  from  that  the 
convection  velocities  Um  and  Uy  of  their  upstream  and  downstream  fronts,  upstream  and 
downstream  being  defined  with  reference  to  the  direction  of  increasing  Z  .  Fig  13b 
illustrates  this  type  of  calculation  for  two  spots  which  were  followed  between  films  5  and  7. 
The  results  are  given  in  Fig  14  for  case  H:  the  ratios  U m/Woo  and  Uy/Woo  are  show  there 
as  a  function  of  Qoo  . 

The  lengthening  of  the  spots  indicates  of  course  that  the  downstream  front  is  moving 
faster  than  the  upstream  front,  but  the  two  velocities  remain  close  at  around  0.7  Woo :  this 
value  is  in  excellent  agreement  with  the  measurements  of  Gaster7,  who  studied  the 
propagation  of  artificial  spots  along  the  leading  edge  of  a  symmetrical  wing.  Beyond 
Qoo  =  70  ms*1  ,  the  spots  become  so  numerous  that  they  combine,  preventing  the 
measurement  of  their  individual  convection  velocity. 

When  considering  studies  in  which  the  flow  is  made  up  of  a  succession  of  turbulent 
spots  and  laminar  zones,  an  important  parameter  is  the  intermittency  factor  y ,  the  fraction 
of  the  total  time  during  which  the  flow  is  turbulent.  The  variation  of  y  as  a  function  of  Z 
was  calculated  in  case  H  for  several  values  of  the  velocity  Q<„  between  64  and  80  ms1  . 
The  results  are  shown  in  Fig  15a,  where  Z  =  0  corresponds  to  the  wing/floor  junction  of 
the  working  section.  At  hot  film  (Z  =  12  cm),  y  remains  equal  to  unity.  After  that,  the 
intermittency  factor  decreases  as  far  as  films  3  and  4  (Z  =  40  to  70  cm),  before  increasing  to 
the  extremity  of  the  wing.  These  results  confirm  the  observations  carried  out  on  the 
fluctuating  signals. 


_ i _ i  - - - — , _ i _ . _ l 

Z(m)  2  0  '  Z(m)  2 


a)  Experiment  b)  Equation  ('5) 

Fig  15  Intermittency  factor  along  the  attachment  line,  case  H 

Stewart  and  Poll15  have  recently  proposed  a  semi -empirical  equation  showing  the  variation 
of  y  on  the  attachment  line  as  a  function  of  the  parameters  R  and  q  : 
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This  formula  was  derived  from  experiments  in  which  contamination  was  induced  by 
cylindrical  wires  surrounding  the  leading  edge  and  attached  perpendicularly  to  the 
attachment  line,  the  distance  Z  being  measured  from  the  wire.  The  theoretical  results  given 
by  (5)  are  shown  in  Fig  15b  for  the  same  conditions  as  Fig  15a.  The  essential  differences 
between  the  two  families  of  curves  is  apparent  at  the  low  values  of  Z ,  where  in  the  current 
experiments,  y  decreases  from  unity;  by  contrast,  in  the  Stewart-Poll  formula  it  increases 
from  a  non- zero  but  very  low  value.  This  is  related  to  the  means  by  which  contamination  is 
brought  about,  ie  a  completely  turbulent  boundary  layer  in  our  measurements,  and  a 
cylindrical  wire  in  the  experiments  which  served  to  establish  equation  (5).  On  the  other 
hand,  in  the  region  in  which  y  increases  (Z  >  0.7  m),  there  is  very  acceptable  agreement 
between  Figs  15a  and  15b,  especially  at  Q«,  ^  70  ms'1 . 

3.5  Second  series  of  hot  Him  measurements 

The  next  item  of  interest  was  the  propagation  of  spots  in  the  X  direction,  with  the 
hot  films  remaining  in  the  same  positions  in  Z  as  before  (Fig  12),  but  displaced  along  the 
upper  and  lower  surfaces  of  the  wing.  The  new  X  positions  are  shown  schematically  in 
Fig  16,  for  case  H;  on  both  surfaces  of  the  wing,  the  two  sets  of  four  films  were  now 
positioned  towards  the  maxima  of  the  external  velocity,  at  10-15%  chord  downstream  of  the 
attachment  line.  The  transient  signals  corresponding  to  this  actual  configuration  are  shown 
in  Fig  17a&b,  for  Q«,  =  61  and  70  ms-1  ,  respectively.  On  the  lower  surface  the  shape  of 
the  signals  differs  little  from  that  already  observed  in  the  first  series  of  experiments;  the 
spots  propagate  from  the  attachment  line  with  little  modification  to  their  structure.  On  the 
upper  surface,  on  the  other  hand,  individual  turbulent  bursts  are  no  longer  seen;  at  61  ms-1 
the  signals  remain  of  a  laminar  nature,  and  at  70  ms-1  fluctuations  of  low  amplitude  become 
apparent. 

These  observations  suggest  that  on  the  upper  surface  the  negative  pressure  gradient 
is  sufficiently  great  to  suppress  the  turbulent  spots  originating  from  the  attachment  line,  at 
least  up  to  a  certain  value  of  Qx ,  and  relaminarisation  may  be  referred  to. 

For  two-dimensional  flow.  Launder  and  Jones16  use  the  parameter  K  ,  to 
characterise  possible  relaminarisation  in  a  negative  pressure  gradient,  defined  by; 


K  = 


dU 

e 

2  dX 


e 


(6) 


For  values  of  K  greater  than  10‘6,  experiment  shows  that  a  turbulent  flow  may 
revert  to  a  laminar  condition.  It  would  seem,  according  to  Beasley17,  that  a  value  greater 
than  about  5  x  10*6  may  be  necessary  for  the  boundary  layer  to  return  effectively  to  this 
state. 
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Lower  surface  Lower  surface 

X' (arbitrary  scale)  Film  No  X' (arbitrary  scale)  Film  No 


0  0.04  0.08  0.12  t(s)  0.16  0  0.04  0.08  0.12  t(s)  0.16 

a)  Q„  =  61ms-1  b)  C^lms-l 

Fig  17  Signals  obtained  from  hot  films,  positions  as  in  Fig  16 

In  our  experience,  this  kind  of  information  does  not  exist  for  three-dimensional 
flow.  One  could  put  forward  the  hypothesis  that  the  critical  values  given  above  remain  valid 
under  the  conditions  of  calculation  of  the  parameter  K  along  the  external  streamline,  ie: 
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-0.2  -0.1  x/c  0  -0.2  -0.1  x/c  0 


Fig  18  Variation  of  the  parameter  K  ,  Fig  19  Position  of  hot  films 
case  H,  Qoc  =  60  ms-1 

Fig  18  shows  the  experimental  distribution  of  K  under  the  conditions  of  case  H, 
for  Qoo  =  60  ms1  (the  velocity  at  which  contamination  commences  on  the  attachment  line). 
On  the  lower  surface,  the  value  of  K  remains  very  small,  while  on  the  upper  surface  it 
reaches  10  x  106,  twice  the  minimum  value  required  for  relaminarisation.  As  K  is 
inversely  proportional  to  Qoo ,  the  criterion  indicates  that  a  velocity  in  the  working  section  of 
at  least  120  ms-1  would  be  required  in  order  not  to  suppress  further  the  leading  edge 
turbulence  in  the  accelerated  region  of  the  aerofoil.  The  detailed  verification  of  the 
relaminarisation  criterion  will  be  the  subject  of  later  experiments. 

This  kind  of  measurement  was  carried  out  only  for  case  H.  However,  calculation  of 
the  range  of  K  showed  that,  at  the  onset  of  contamination  its  maximum  value  lies  between 
8  and  10  x  IQ-*6  for  "high  incidence"  configurations  (cases  D  to  H),  and  that  it  is  in  the 
region  of  6  x  lO-6  for  "low  incidence"  configurations  (cases  A,  B  and  C).  Relaminarisation 
would  therefore  seem  possible  for  all  the  combinations  of  a  and  <p  considered  in  this 
study.  The  high  values  of  K  are  attributable  to  the  shape  of  the  leading  edge  which  induces 
particular  velocity  distributions;  it  is  evident  in  Fig  1  la&b  that  the  concavity  of  Ue(X)  faces 
upwards,  it  on  the  upper  surface  there  are  velocity  gradients  more  intense  than  in  the  region 
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of  the  attachment  line.  With  "classical"  aerofoils,  on  the  other  hand,  the  greatest 
accelerations  are  measured  on  the  attachment  line  and  relaminarisation  cannot  occur 
downstream  at  Reynolds  numbers  at  which  contamination  appears,  as  demonstrated  by 
Poll's  analysis  for  elliptic  leading  edges9. 


Lower  surface 

X'  (arbitrary  scale)  Film  No 


I 


Fig  20 


Upper  surface 


Signals  obtained  from  the  hot 
films,  positions  as  In  Fig  19, 
Q„  =  64  ms*1 


The  hot  films  were  finally  removed  and  redistributed  along  the  chord  between  the 
attachment  line  and  13%  chord;  their  new  arrangement  is  shown  schematically  in  Fig  19, 
while  the  recordings  obtained  at  Q*,  =  64  ms*1  are  shown  in  Fig  20.  Two  forms  of 
transition,  totally  different  in  nature,  are  observed;  the  lower  surface  is  swept  by  turbulent 
spots  originating  in  the  attachment  line.  These  are  suppressed  on  the  upper  surface  (film  4) 
where  the  transition  occurs  downstream  in  the  region  of  positive  pressure  gradient,  by  a 
streamwise  instability  (films  6  and  8). 


4  CONCLUSION 

This  study  has  allowed  a  certain  amount  of  information  to  be  gathered  on  the 
mechanisms  of  transition  in  three-dimensional  incompressible  flow. 

The  prediction  of  transition  induced  by  crossflow  instability  is  relatively  easy  in 
regions  of  intense  negative  pressure  gradient,  for  example,  between  the  leading  edge  and  the 
point  of  maximum  velocity.  Such  a  prediction  may  be  made  by  the  en  method,  as  well  as 
by  the  use  of  simple  empirical  criteria.  In  the  case  of  a  wing  having  a  cambered  leading 
edge,  as  considered  in  the  first  series  of  experiments,  the  problems  become  complicated 
downstream  of  the  maximum  velocity  where  pressure  gradients  remain  low.  The 
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appearance  of  crossflow  velocity  profiles  in  S  and  the  simultaneous  development  of 
streamwise  instability  make  it  necessary  to  fall  back  on  the  theory  of  laminar  instability. 

Regarding  the  problem  of  leading  edge  contamination,  it  has  been  shown  in  the  two 
series  of  measurements  that  the  critical  value  R  =  245  may  be  quite  confidently  employed  to 
predict  the  Reynolds  number  beyond  which  turbulent  spots  may  develop  along  the 
attachment  line.  However,  before  this  critical  value  is  achieved,  a  significant  part  of  the 
surface  of  the  model  may  be  swept  by  turbulent  bursts  originating  from  the  junction  of  the 
wing  and  fuselage.  It  must  also  be  borne  in  mind  that  contamination  phenomena  are  not 
uniform  along  the  span  (y  varies  in  the  Z  direction)  and,  on  account  of  this,  a  calculation 
of  the  boundary  layer  could  not  be  carried  out  using  the  classical  hypothesis  dfdZ  =  0 ,  even 
if  the  distribution  of  external  velocity  does  not  vary  much  in  the  direction  parallel  to  the 
leading  edge. 

These  experiments  have  also  shown  that  relaminarisation  would  appear  to  be 
possible  if  the  pressure  distribution  downstream  of  the  attachment  line  possesses  streamwise 
gradients  which  are  greater  than  along  the  attachment  line  itself. 
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